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We find that sudden future singularities of pressure may also appear in spatially inhomogeneous 
Stephani models of the universe. They are temporal pressure singularities and may appear indepen- 
dently of the spatial finite density singularities already known to exist in these models. It is shown 
that the main advantage of the homogeneous sudden future singularities which is the fulfillment of 
the strong and weak energy conditions may not be the case in the inhomogeneous case. 
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I. INTRODUCTION 

It has been shown by Barrow 1] that for Friedmann 
cosmological models which do not admit an equation of 
state which links the energy density g and the pressure 
p, a sudden future (SF) singularity of pressure may ap- 
pear, even for the matter fulfilling the strong energy con- 
dition g > 0, g + 3p > 0, though violatingthe domi- 
nant energy condition g > 0,—g < p < g This is 
in contrast to the most of the current observational dis- 
cussion in cosmology, mainly devoted to determining of 
the barotropic index w in a barotropic equation of state 
p — wg, which tightly constraints the energy density and 
the pressure (3(. On the other hand, the observational 
data interpreted by such an equation of state cannot ex- 
clude a possibility of barotropic phantom cosmological 
models :4|. These models violate null energy condition 
g+p > 0, and consequently all the remaining energy con- 
ditions 5|. Besides, phantom models allow for a Big- Rip 
(BR) curvature singularity, which appears as a result of 
having the infinite values of the scale factor a(t) at finite 
future. This is in opposition to a curvature Big-Bang 
(BB) singularity which takes place in the limit a — > 0. 

The common feature of BB and BR singularities is 
that both g and p blow up equally. This is not the case 
with a SF singularity for which a blow up occurs only 
for the pressure p, but not for the energy density g. It is 
interesting that SF singularities are similar to those ap- 
pearing in spatially inhomogeneous Stephani models of 
the universe [|J, in which they were termed finite den- 
sity (FD) singularities @>@- However, FD singularities 
occur as singularities in spatial coordinates rather than 
in time (as SF singularities do), which means that even 
at the present moment of the evolution they may exist 
somewhere in the Universe 0, 0, . In this letter we 
show that sudden future (SF) singularities (as temporal 
singularities) can be inhomogenized in the sense, that 
they may appear in spatially inhomogeneous models of 
the universe, independently of the spatial finite density 
(FD) singularities allowed in these models. We also show 
that the inhomogeneous Stephani models lead to energy 



conditions violation, which mainly refers to the fact that 
they admit FD singularities. 

The SF singularities appear in the simple framework 
of Friedmann cosmology with the assumption that the 
energy-momentum is conserved, so that one can write 
the energy density and pressure as follows (following pj 
we assume 8ttG = c = 1, k = 0, ±1) 
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From J2I one is able to notice that the singularity of 
pressure p — > oo occurs, when the acceleration d — > — oo. 
This can be achieved for the scale factor 



„(/) = . i + („„- i.(ij -Afi-i 



(3) 



where a s = a(t s ) with t s being the SF singularity time 
and A,q,n — const. It is obvious from that <z(0) = 
and so at zero of time a BB singularity develops. For the 
sake of further considerations it is useful to write down 
the derivatives of the scale factor ©, i.e., 



a = qt s {a a - A) 



= q(q-l)ti(a s -A) ^- 



, n q— 1 / , \ n— 1 

q-2 



, (4) 



(5) 



The main point is that the evolution of the Universe, as 
described by the scale factor |J5J, begins with the stan- 
dard BB singularity at t = 0, and finishes at SF singu- 
larity at t — t s , provided we choose 



1< n < 2, i)<q<l 



(6) 
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For these values of n and q, the scale factor (0 vanishes, 
and its derivatives ©-© diverge at t — 0, leading to a 
divergence of g and p in (BB singularity). On the 

other hand, the scale factor © and its first derivative (QJ 
remain constant, while its second derivative JSJ diverge, 
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leading to a divergence of pressure in (J2J only, with finite 
energy density fl). This behaviour means, for example, 
that positive curvature (fc = +1) Friedmann models may 
not recollapse to a second BB singularity - instead they 
terminate in a SF singularity [l2|. 



II. INHOMOGENEIZED SUDDEN FUTURE 
SINGULARITIES 

Now, let us consider inhomogeneous Stephani mod- 
els. They appear as the only conformally flat perfect- 
fluid models which can be embedded in a 5-dimensional 
flat space 0|. Their metric in the spherically sym- 
metric case reads as (notice that we have introduced a 
Friedmann-like time coordinate which eliminated one of 
the functions of time in the original Stephani metric §) 
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where 



V(t,r) = 1 + -Mt)r 2 



')] , (7) 



(8) 



and (...)' = d/dt. The function a(t) plays the role 
of a generalized scale factor, k(t) has the meaning of a 
time-dependent " curvature index" , and r is the radial 
coordinate. 

Their analogy to SF singularity models is that they do 
not admit any equation of state linking p to g through- 
out the whole evolution of the universe, although at any 
given moment of the evolution, such an equation of state 
(varying from one spacelike hypersurface to the other) 
can be admitted. An analytic equation of state can also 
be admitted at any fixed subspace with constant radial 
coordinate r, but not globally Q. For the sake of simplic- 
ity, first the spherically symmetric models will be consid- 
ered (note that in [H, H a time coordinate t analogous to 
the cosmic time in Friedmann models was marked by r, 
which had nothing to do with a common conformal time 
coordinate). The energy density and pressure are given 
byi 



g(t) = 3C 2 (t) 



r d 2 (£) Jfe(t) 



a 2 (t) 



2 (t) 



p(t,r) = - 3C 2 (t) + 2C(t)C(t) 
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(9) 
(10) 



and generalize the relations (JTJ and J5J to inhomogeneous 
models. 

We now show that it is possible to extend SF singulari- 
ties into inhomogeneous models. Following p| we assume 
that the functions k(t) and a(t) are related by 



k(i) = -aa(t) 



(11) 



with a = const. In fact, the limit a — > gives the Fried- 
mann models (cf. the discussion of the conditions to 
derive such a limit in |8|]). Inserting into IjlUI) we get 
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From Ijl-'SH one can see, that p — > oo, when acceleration 
d — * — oo for an arbitrary value of the radial coordinate 
r. We can then say that we generalized SF singulari- 
ties (given, for example, by the scale factor 0) onto an 
inhomogeneous model of the universe. 

However, it is very interesting to notice that in such 
a generalization not only SF singularities appear, but 
also FD singularities are possible for the radial coordinate 
r 2 — > oo. These seem to be far away from us, and so not 
very harmful, since r 2 — > oo defines an antipodal center of 
symmetry in the spherically symmetric Stephani models. 
The advantage of these FD singularities is that they are 
able to drive the current acceleration of the universe [Tol 

mm 

The procedure of inhomogenizing SF singularities may 
be extended into the general Stephani models for which 
there is no spacetime symmetry at all, and so they are 
completely inhomogeneous. The most general Stephani 
metric in cartesian coordinates (x. v. z) |3.[l^ reads as 



ds 2 
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(14) 



where 



V(t,x,y,z)= (15) 
1 + h(t){[x - x (t)} 2 + [y- y {t)} 2 + [z- z {t)} 2 } , 

and Xo,yo,zo are arbitrary functions of time. Now the 
general expression for the pressure is (the expression for 
the energy density © remains the same) 



p(t,x,y,z) = - 3C 2 (t) + 2C(t)C(t) 
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Inserting the time derivative of 10 and the function 
V(t, x, y, z) from ltT5|l into l(TB|) gives 
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(17) 
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It is easy to notice that SF singularity p — > ±oo appears 
with © for d — > — oo, if (V/a)/(V/a)' is regular and the 
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sign of the pressure depends on the signs of both a/ a and 
(y/a)/(V/a)' . This proves that we can inhomogenize SF 
singularities for a Stephani model with no symmetry. 

In fact, SF singularities appear independently of FD 
singularities whenever a — > — oo and the blow-up of p is 
guaranteed by the involvement of the time derivative of 
the function C(t) in ifTUI) . 

That makes a complimentary generalization to the one 
which extends SF singularities into the theories with 
actions being arbitrary analytic functions of the Ricci 
scalar and into anisotropic (but homogeneous) models 

It appears that the main motivation for studying SF 
singularities Q was the fact that, unlike phantom mod- 
els |4|, they obey the strong and weak energy conditions, 
though they do not obey the dominant energy condition. 
In this paper we raise the point that the question of pos- 
sible violation of the energy conditions for the inhomog- 
enized SF singularity models is a bit more complex than 
for the homogeneous ones. From (|1L)|) and i|13|) for the 
strong, weak and dominant energy conditions to be ful- 
filled we have: 
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(18) 



(19) 



(20) 



In fact, the dominant energy condition requires fulfilling 
both l(T§|) and (|2*U|) . Notice that in view of © the energy 
density in inhomogeneous Stephani models is always pos- 
itive, i.e., 



g> 



(21) 



This means that the strong and weak energy conditions 
are still not violated if d — > — oo, in analogy to homoge- 
neous models, provided 



1 a 9 
a 4 



(22) 



since the first term with a in l|18|) and (|19(l must dominate 
the second (remember that d — > — oo). Notice that the 
equality 1/a — ar 2 /4 may lead to a pressure singularity 
avoidance in l|13l) . Assuming that the generalized scale 
factor a(t) > 0, we conclude from lll'L'1) that the strong 
and weak energy conditions are always fulfilled, if a < 0. 
However, an accelerated expansion for an observer at the 
center of symmetry at r = can only be achieved, if 



a > 18] (the spatial acceleration scalar reads as ii — 
— 2ar - lower pressure regions are away from the center). 
This means that the strong and weak energy conditions 
are not necessarily fulfilled for the models with a > 0. 
In particular, they cannot be fulfilled at an antipodal 
center of symmetry at r 2 — > oo, unless a — > (where Big- 
Bang singularity appears and so BB and FD singularities 
coincide - see Section ITTT|) . 

On the other hand, the first part of the dominant en- 
ergy condition may not be violated if the contribution 
from the last term with d in l|2U|) . which includes r 2 , 
does not overweigh the first one, i.e., when 
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(23) 



This should be appended by the condition i|19fl which is 
equivalent to l|22|l . i.e., the dominant energy condition is 
fulfilled if 



1 a 9 1 
- < —r z < - 

a 4 



a 



(24) 



This is obviously a contradiction which means that, sim- 
ilarly as in the isotropic Friedmann models, SF singular- 
ities violate the dominant energy condition. 

Such a violation of the dominant energy condition 
also appears in M-theory-motivated ekpyrotic models in 
which p 3> g during recollapse [l9| . 

Let us now discuss the problem of the possible energy 
conditions violation in the general Stephani model. Using 
© and (|17|) we get for the strong, weak, and dominant 
energy conditions 
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Obviously, the dominant energy condition requires fulfill- 
ing both (|26|l and 1)27(1. Before we go any further, using 
IjlSp. we note that 



V(t,x,y,z) 



a(t) 
1 k 
'la 



(28) 



(x - x ) 2 + (y- y ) 2 + (z- z ) 2 
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V(t,x,y,z) 



a(t) 
a 1 
a 2 4a 



k-k- 



(29) 



(x - x Q ) 2 + (y- y ) 2 + {z- z Q ) 2 



k 

[(x - x ) x + (y - yo) Vo + {z- z ) z ] . 

It is important to notice that the ratio of l|28|) and l|29|) 
which appears in the conditions Q25J1. I|26|) . and l|27|) al- 
lows to cancel 1/a from both the numerator and the de- 
nominator. Apart from that, one is able to take a/a out 
in (|29() and cancel it with the same term standing in front 
of the last term in these conditions. This basically means 
that, bearing in mind the fact that a — > — oo, the strong 
and weak energy conditions are fulfilled provided one of 
the expressions 
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ka 



- [{x - x ) x + (y- yo) Vo + (z- z ) z ] - 1,(31) 
a 

is negative. It is clear that for k(t) > (|30|l is always 
positive, so that l|31|) must necessarily be negative and it 
certainly does, at least in some regions of space. On the 
other hand, for k(t) < <|3l)|l can be both positive and 
negative which requires l|31[) to be negative and positive, 
respectively. In conclusion, similarly as in the spheri- 
cally symmetric case, the strong and weak energy condi- 
tions may be violated for inhomogenized SF singularities 
in Stephani models. This is different from what we have 
for isotropic SF singularities. Finally, one can easily no- 
tice that in order to fulfill the dominant energy condition 
the ratio of (|30|) to l|31|) should be simultaneously positive 
and negative, which is a contradiction. This means that, 
like in the isotropic case, a general Stephani model al- 
lows SF singularities which always violate the dominant 
energy condition. In conclusion, one can say that the 
problem of the energy conditions violation by SF singu- 
larities in the inhomogeneous models is more complicated 
than in the isotropic ones and so the results based only 
on the isotropic models cannot be trusted. 



III. INHOMOGENEOUS FINITE DENSITY 
SINGULARITIES 



Under the choice of l|ll(l . we have (k/a)' = 0, and so the 
singularities appear at r 2 — > oo. In general, it may not 
be so, which was explicitly shown in y|. For example, by 
choosing 



k(t) 



a(t) = at 2 + /3< + 7 
- a 2 = -4aa(t) + A 



with 



A = 4a 7 + 1 - /3 2 , 
the FD singularities appear for @ 

|r|=2/V=A. 



(33) 
(34) 



(35) 
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Of course the condition 111|) is obtained in the limit 
A — * from l|34|) which moves FD singularities to an 
(30) antipodal center of symmetry at r 2 — * oo. Having cho- 



sen 7 = 0, A = 1 - 1 in (Model I of Ref. [ljj, 
called Dabrowski model in Ref. [ll|) the energy density 
and pressure are then given by 
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For the strong, weak and dominant energy conditions to 
be fulfilled, respectively, we get the requirements 



q + 3p = —6a 



t(at + (3) 



>0 , 



V = 2 l -^ at+ l )r2 >0 



Q-p 
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(39) 
(40) 
(41) 



If a > (acceleration then the strong energy con- 

dition is violated if 



t(at + 0) > 



(42) 



and this may happen independently of the radial coor- 
dinate r. The weak energy condition is violated for the 
domain of space in which 



~2 > at(at + [3) . 



(43) 



In the context of temporal SF singularities of pressure 
we will further briefly discuss the occurrence of spatial 
FD singularities of pressure in the Stephani models and 
possible energy conditions violation. From IjlQI) we can 
see that FD singularities appear whenever the radial co- 
ordinate 



= -4 



(a) 



(32) 



With the strong energy condition violated, this gives a 
weak energy condition violation only in some spatial do- 
main since the right-hand side of 1)43(1 has a positive value, 
but including the center of symmetry at r = 0. On the 
other hand, for the decelerated expansion, a < 0, and the 
right-hand side of (|43() has a negative value, so the weak 
energy condition is violated everywhere in the universe 
(i.e., for all values of r). If the strong energy condition is 
not violated and a > (acceleration), then again there 
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is a weak energy condition violation for all values of r. If 
the strong energy condition is violated, and a < 0, then 
the weak energy condition is violated only in some spatial 
domain which, however, includes the center of symmetry 
at r = 0. The dominant energy condition is violated for 

-1 < -lat(at + 0) , (44) 
r z 2 

which, combined with l|43(l . gives 

at{at + (3) < 4r < -\at(at + /3) . (45) 
r 2 - 2 

This last condition can only be fulfilled either if a < 
and t(at + (3) > 0, or if a > and t(at + /?) < 0. Then, 
at least for these particular class of Stephani models, FD 
singularities may lead to a violation of the energy con- 
ditions in a similar way as BR singularities in phantom 
cosmology do. 

An interesting problem is a possible avoidance of SF 
and FD singularities in the universe. SF singularities can 
easily be avoided by imposing an analytic form of the 
equation of state p = p(g). Even without this assump- 
tion, some other ways of their avoidance by introduc- 
ing quadratic in Ricci curvature scalar terms 0, or by 
quantum effects |2(]|. are possible. On the other hand, a 
necessary condition to avoid FD singularities in Stephani 
models comes from ijTUj) and reads as 

jty > • ( 46 ) 

In our special model they can be simply 



avoided, if 

A > . (47) 



IV. CONCLUSIONS 

In conclusion, we have shown that one is able to spa- 
tially inhomogenize sudden future (SF) singularities in 
the sense that these singularities do appear in inhomo- 
geneous models of the universe. However, despite ho- 
mogeneous SF singularities, they may violate the strong 
and weak energy conditions in some regions of space, al- 
though they share the dominant energy condition viola- 
tion with homogeneous models. It shows that making any 
important conclusions about physics, on the basis of the 
isotropic models only, may be misleading and should not 
be trusted. A possible violation of all the energy condi- 
tions by inhomogenized SF singularities is similar to what 
happens to Big- Rip singularities in phantom cosmologies. 

Besides, we have noticed that, apart from sudden fu- 
ture singularities, the inhomogenized models also admit 
finite density (FD) singularities which are spatial rather 
than temporal. In relation to this we have discussed an 
example of an inhomogeneous model with spatial finite 
density singularities of pressure and studied the domains 
of its energy conditions violation. 
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